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Abstract. For the hyperbolic system of quasilinear first-order partial differential equations, 
linearizable by hodograph transformation, the conservation laws are used to solve the Cauchy 
problem. The equivalence of the initial problem for quasilinear system and the problem for 
conservation laws system permits to construct the characteristic lines in domains, where 
Jacobian of hodograph transformations is equal to zero. Moreover, the conservation laws 
give all solutions of the linearized system. Some examples from the gas dynamics and theory 
of plasticity are considered. 
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1 Introduction 

The Riemann method is one of the classical ones and is used widely for solution of Cauchy 
problem for hyperbolic linear equation of second order for one function of two independent 
variables 

d*^ u du du 

+ a{x,y)— + b{x,y)— + c{x,y)u = f{x,y), (1) 



dxdy ^ dx ' 9y 

where u{x,y) is specified on the boundary L = {{x,y)\x = x{T),y = ?/(t)}, which is not 
a characteristic line. 

To determine the value of u at point x = x' , y = y' , the so-called Riemann function 
R{x,y; x' ,y') is used, which is a solution of the special characteristic Cauchy problem (Goursat 
problem) for the homogeneous adjoint equation, corresponding to (1) 

#B 3(aR) Om^^^^^^ 



dxdy dx dy 



0' ^\x=x',y=y' — 1- 



y=y 



If the Riemann function is determined, then the solution of Cauchy problem for (1) is given by 
an integral representation. The main reason for the introduction of adjoint equation (2) is to 
make the line integral vanish around closed paths [15]. In other words there is a conservation 
law of the special form, that, for example in [6] is called the obvious one. 



*This paper is a contribution to the Special Issue "Geometrical Methods in Mathematical Physics". The full 
collection is available at http://www.emis.de/journals/SIGMA/GMMP2012.html 
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The same method is apphed for linear system obtained from quasilinear one by hodograph 
transformation. This 'speedgraph' transformation, see for example [28], is just an interchange 
of roles of the unknown functions and the independent variables. There is a lot of works, where 
the hodograph-type transformations are used to linearize different kinds of systems of different 
orders. Let us mention paper [7], where an algorithm for establishing whether a given quasilinear 
equation is linearizable (is solvable in terms of either a linear partial differential equation or of 
a linear integral equation) is proposed to apply Painleve test to quasilinear equations. The 
so-called extended hodograph transformation (including integral term) is introduced to reduce 
quasilinear PDE to the semi- linear one. In work [9] the hodograph transformation with the 
principle of nonlinear superposition is used to produce some new solutions of the special classes 
of equations. 

But the linearization of quasilinear systems considered later by hodograph transformation is 
valid only in the domains, where the Jacobian of transformation is not equal to zero, which is un- 
known before the solution of the quasilinear system, because of the dependence of characteristic 
curves on the solution. 

Conservation laws of quasilinear homogeneous hyperbolic system for two functions of two 
independent variables, related to the solution of corresponding linearized system are used in the 
paper for the solution of the Cauchy problem. Vanishing of Jacobian is not a restriction now, 
this allows to construct the characteristic fields corresponding to the simple waves. 

The paper is structured as follows. In Section 2 we stand the problem and describe some 
basic properties of a hyperbolic quasilinear system. Section 3 deals with the construction of the 
conservation laws of the special form, its relation with a solution of the linearized system and 
the description of the solution of Cauchy problem. Some applications of the exposed method 
are considered in Section 4, in particular the Cauchy problem from the theory of plane plasticity 
for the loaded cavity is solved for any convex form of contour. 



2 Hyperbolic quasilinear system 

Let us consider a quasilinear system of homogeneous PDEs of two independent variables x, y 
and two dependent ones u, v in the form [21] 

dx ~^ dy ^ ^ ^ 

where A = \\aij{u,v)\\, B = \\bij{u, v)\\, i, j = 1,2, U = {u, v)"^. 

If matrix A is not degenerate, then system (3) can be written in the normal form 

dU ^dU _ ^ 
dx dy ' 

where M = \\mij{u,v)\\. 

Let us set up the Cauchy problem for system (4) due to [21]: in some neighborhood of the 
arc C: a ^ r ^ 6 of an initial curve L = {{x,y) : x = x{T),y = ?/(t)} it is necessary to determine 
the solution of (4) satisfying the initial condition on L 

U{x{T),y{r)) = U\T), rG[a,6]. (5) 

The Cauchy problem (4), (5) is supposed to be the normal one (in the sense of [21]) and the 
solution exists in some neighborhood of L. 

Let system (4) be a strictly hyperbolic one. It means that matrix M has two real different 
eigenvalues Ai and A2, obtained as roots of the equation 

det(M - \E) = =^> 2Ai^2 = m-ii + "^22 ± \/ ("in - 77122)^ + 4mi2m2i, 

that gives two left eigenvectors l\ = and I2 = (^25^2) respectively. 
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Let us consider the differential forms 

(jjfc = ll{u, v)du + l\{u, v)dv = 0, = 1, 2, 

which can be integrated, because in this case there always exist integrating factors. Then, the 
corresponding two integrals ^^{u^v) = const can be taken as Riemann invariants = ^k{u,v) 
and system (4) takes a diagonal form 

dR dR 

— +A— = 0, A = diag(Ai,A2), R={ri,r2f. (6) 
System (6) has two families of real characteristic curves determined by the following equations 

^ = Ai(n,f), ^ = A2(u,i;), Ai / A2, (7) 

that is, variable ri is the invariant along the first characteristic curve and r2 is the invariant along 
the second one. For more detailed theory of Riemann invariants see, for example, [12, 20, 21]. 

The systems of such a form are widely used in the mechanics of a continuum media [21], 
in the gas dynamics for describing isoentropic plane-symmetry flows; in the theory of plane 
plasticity for the stresses of a deformed region under the different yield criterions [13], for the 
motion of granular materials, for the propagation of the plane wave of loading in homogeneous 
semi-infinite elastic-plastic beam [18], etc. 

It is well known, that system (4) can be linearized by a so-called hodograph transformation 
of the form x = x{u, v), y = y{u, v). Thus, system (3) takes a linear form 

dei{A){A^)~\uy{u,v) = det{B){B^)~^Vux{u,v), Vu = id/du,d/dvf , 

and for diagonal form (6) one can obtain 

AV/jy(ri, r2) = det{A)V Rx{ri, r 2), Vr = {d/dri, d/dr2f . (8) 

Let us note, that it is possible to obtain the solution of system (6) from the solution 
of (8) and vice versa only when two corresponding Jacobians Ji = \d{ri,r2)/d{x,y)\ and 
J2 = \d{x,y)/d{ri,r2)\ are not equal to zero 

J dri dr2 dr2 dri 

dy dx dy dx ^ 

Let us mention, that all solutions of system (6) can be divided into two classes: the singular 
and non-singular ones. The singular solutions are obtained as solutions of equations Ji = 0. 
And non-singular ones are the solutions of system (8). Solutions with degenerate Jacobi matrix 
form a class of solutions called multiple waves and can be found by so-called method of the 
degenerate hodograph (see, for example [16]), and in general belong to the class of partially- 
invariant solutions [19]. A solution of considered system (6) for which the rank of the Jacobi 
matrix in a domain G G satisfies the condition 

d{ri,r2) 
rank ■ 



d{x,y) 



is called a multiple wave of the rank r. In general r < min {n,m}, where n is a number of 
independent variables and m is a number of functions. There are only two possibilities for the 
rank in our case: r = Oor = l (r = 2 corresponds to non-singular solution). The case r = 0, 
means that both Riemann invariants are constant functions, that corresponds to the constant 
functions u and v. If r = 1, then a multiple wave is called a simple wave and is widely known. 
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Figure 1. Areas of straight line cliaracteristics. 



For this kind of solution there are finite relations between the functions u, v. For example, 
simple waves for the plane plasticity equations, considered later, were analyzed in [14]. 

The main difficult in practical applications is the determination of a domain G, mentioned 
above, because the domain of determinacy of the Cauchy problem for the system of quasilin- 
ear equations is determined simultaneously with the solution and, generally speaking, cannot 
be indicated beforehand [21]. In other words, before the complete solution of Cauchy prob- 
lem (4), (5) it is impossible to calculate both Jj. By the same reason, before the complete 
solution it is difficult to point out the domain, where the characteristic curves (7) do convert to 
straight lines (see Fig. 1) for given form of the initial curve and for given initial conditions. 

Due to the linearity of (8) any non-singular solutions can be obtained from the one by the 
superposition principle (that corresponds to the symmetry, admitted by a linear system) . In the 
theory of group analysis such kind of systems is called automorphic: if one particular solution 
is known, then the orbit of this solution under the group of admitted symmetries forms locally 
a general solution of the given system. In such a way, for example, some new exact solutions 
for the plane plasticity system were constructed in [29] . Using conservation laws permits to find 
out both nonsingular and singular solutions. 

System (4) can be extended for three or more functions and in this case is called one- 
dimensional system of hydrodynamic type [27]. A natural hamiltonian formalism was proposed 
for this class of homogeneous systems of PDF and the generalized hodograph method generates 
from every its solution a symmetry (commuting flow) that finally leads to solution. Another 
approach to generalize the concept of Riemann invariants to systems with more than two inde- 
pendent variables (m > 2) or with the number of functions more than two (n > 2), is related 
with so-called simple integral elements [10]. 

3 Conservation laws 

The concept of conservation laws is one of the fundamental characteristics of a real physical 
process. One can find the formal definition of the conservation laws of the system of diffe- 
rential equations, for example in [3], where the concept of the so-called operator of universal 
linearization is used. 

Let us seek the conservation law of system (6) directly in the form 




(9) 



which is valid for any solution of (6). Let us multiply (6) by a vector a = (01,02) [21], then. 
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Figure 2. Cauchy problem. 



eliminating q, we obtain the linear system for functions and '0 

AVHV3(ri, ra) = VRV(n, r-a). (10) 

Now let us describe the way to solve Cauchy problem using conservation laws. Let P{x{a),y{a)), 
Q{x(b),y{b)) be two end-points of the arc C, M{Mx,My) be a point of intersection of two 
characteristic lines: ri = r^, going from the point Q and r2 = rg, going from the point P (see 
Fig. 2). 

Let us apply the Riemann method not for linearized system (8), but for system (10). Taking 
integral over the closed path PQM, which due to Green theorem is equal to zero, we have 



il^dx — (fdy = = / tpdx — (fdy + / (il; — (fXi) dx + / (■0 — (^Aa) dx. 

PQM JPQ Jri=rO Jr2=r0 

Integrating by parts the above two integrals 

/ {ip - ifXi) dx = {it - ifXi) dx = x {ip - ipXi) HZq'' - xdi^p - ifXi) , 

Jri=r° JqM JqM 

{ip-ipX2)dx= / {tp - ipX2)dx = xiip - (pX2)\IZm ~ / xd{ip - (pX2) , 

r2=rq J MP "^-^^^ J MP 

and taking (without loss of generality) the following conditions 

(^ - >^lf)\r,=rO = 1, (V' - ^2^)\r,=r0 = 0, (11) 

we get coordinate 

M^^m-IJ/^-/i)^r. (12, 
Analogically, for the y-coordinate of point M we obtain 

My = y{b) -I f - dr, (13) 



PQ 



dr dr 
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but now the conditions for the functions (p, ij) are the following ones 

(^/Ai - (/.)|,^^,o = 1, (V/A2 - <^)|,,=,o = 0. (14) 

In such a way, the coordinates of point M are known and one can reconstruct the values of 
functions u, v from the initial condition (5). 
Finally, we have demonstrated the following 

Theorem 1. The solution of Cauchy problem for quasilinear system (4), (5) is equivalent to 
the solution of two Cauchy problems for the linear system (10) with conditions (11) and (14). 

The main application of Theorem 1 is that one can obtain the solution of Cauchy problem 
for the original system without any restrictions on the Jacobian of hodograph transformation. 

Note, that the above procedure can be used to solve the characteristic problem, when the 
initial data are prescribed on two intersecting characteristics (see [23]). 

Let us note, that in [6] the classification of second-order linear differential equations with 
two independent variables in terms of first-order conservation laws is proposed. System (6) can 
be reduced to such an equation, and in this case conservation laws (9) are among the first-order 
ones. 

Let us express (8) and (10) in terms of the hyperbolic equations of the second order 

d'^x 1 8X2 dx ^ 1 dXi dx _ ^ ^^^^ 

dridr2 Xi — X2 dr2 dri Ai — A2 dri dr2 

^ 1 ^Ai dip 1 8X2 /-,^g^ 

dridr2 Ai — A2 dr2 dri Ai — A2 dri dr2 

The most general equivalence transformations, preserving the differential structure of these 
equations have the form [19] 

Ri = fM), u{ri,r2) = w{ri,r2)V{Ri,R2), i = 1,2. (17) 

Using this concept, we can demonstrate the following 

Proposition 1. The conservation laws (10) give all nonsingular solutions of linearized sys- 
tem (8). 

Proof. Let us consider the well-known lemma about Laplace invariants (/i^*), /c*^*)) for two equiv- 
alent hyperbolic equations [19], namely, two hyperbolic equations of the form 

^ + C«(n,.2)^ + C«(n,.2)^ + C«(n,r2)y(^) = 0, 

are equivalent iff there exist a functions /j, w{ri,r2) (17), such that its corresponding Laplace 
invariants 



dri ' ^ ^ ' dr2 
dKi dR,2 



are related by the following equalities 

h^'\ri,r2) = f^^h('\fi{ri)j2{r2)), fe«(ri,r2) = ^§^k('\fi{ri)j2{r2)). 
ari ar2 dri dr2 
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In particular, two equations are equivalent up to the factor-function w{ri,r2) only, iff 
/j{i)(ri,r2) = /i(2)(i?i,i?2), fc(i)(ri,r2) = k^^\Ri,R2). 
Comparing the Laplace invariants of (15) and (16) 

hi-) = 1 ^'^2 _ 1 ^^2 dX2 ^ ^(^)^ 

Ai - A2 dridr2 (Ai - A2)^ dri dr2 
^ 1 ^'^1 _ ^^_9Ai aAi ^ ^(^)^ 
Ai - A2 dridr2 (Ai - A2)^ dri dr2 

one can obtain the relation between eigenvalues 

(Ai-A2)^^(Ai + A2) = ——-— — . (18) 
dridr2 dri dr2 dri dr2 ' 

Function w can be determined from the following system of equations 

1 dw _ (^) (^) 1 dw _ (^) (^) 

~Q —^2 ~ ^2 ' ~ "-^l "'-'1 ' 

w ori w or2 

/^ + c{-^'^ + ci-^^ = 0, (19) 
orior2 or I or 2 

where C^''\ cf^ are the corresponding coefficients in equations (15) and (16). 

Let us note, that a similar comparison of Laplace invariants one can make for the equations, 
obtained from (10) and (8), but expressing functions -0 and y. ■ 

In particular, two simplest cases can be picked out: 

1) if det(A) = K = const, then 99 = y{ri,r2), tp = Kx{ri, r2); 

2) if Ai = -A2, then = x{ri,r2), ip = -y(ri, r2). 

An example of the first case is the system of the ideal plane plasticity [13] 



M 



— cot2f — 2A;/sin2w 
-l/(2/csin2u) -cot2u 



with det(A) = -1. 

One of the examples of system (4) for the second case mentioned above, is system 
mil = 77122 = 0, mi2 = -l, m2i=F^(ii), 
related to the well known equation 
d^Y ^2fdY\ d'^Y 



dx^ \dy ) dy"^ ' ^ 

investigated in [30] in the context with Fermi, Pasta and Ulam (1955) results on the vibration 
of a nonlinear, loaded (or beaded) finite string. In such a case we have Ai = — A2 = —F. 
We consider these and some other examples in the next section. 

4 Examples of application 

It is well known, that linear homogeneous hyperbolic equations (1) are completely classified [19] 
with respect to the group of admitted point symmetries. There are three different kinds of such 
equations, depending on the form of its Laplace invariants. As for Goursat problem (2) for 
Riemann function, it can be reduced to the solution of an ordinary differential equation in the 
same classified cases [11]. Let us consider some examples of these kinds of equations. 
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4.1 Plane ideal plasticity with Saint-Venant— Mises yield criterion 

This system was investigated, using the group of admitted symmetries: for its invariant solutions 
see [1], all its conservation laws and highest symmetries were described in [22] and for the 
reproduction of solutions by point transformations see [24, 25, 29]. Being semi-inverse method, 
group analysis provides analytical solutions and then one can determine the boundary conditions 
for obtained solutions. But if the the boundary conditions are given from the beginning, then 
the method of conservation laws can be applied to solve the boundary problem directly. 
The system has the form 

d(j ( do do \ Oct ( QO dO \ 

— -2fc — cos 26* + — sin 26* =0, — - 2A; — sin 26* - — cos 26* = 0, (21) 

ox \ox oy J oy \ox oy J 

where a is a hydrostatic pressure, is an angle between the first main direction of a stress tensor 
and the Ox-axis, A; is a constant of plasticity. 

The functions for this system in the form (6) are as follows 

u = a, V = 0; Ai=tanf, A2 = -cot7;; ^ 2A; ~ ''^^ ^ 2fc ^' 

In this case det(A) = —1, so f = y{ri,r2), ip = —x{ri,r2). 
The solution of problem (10), (11) has a form [23] 

dp , dp 

if = 2- — cosf — psmw, yj = 2— — svuv + pcosv, 
ori ori 

where function p(ri,r2) looks like this 

/5(ri,r2) = (ri,r?,r2,r2) cos l^^-^-^J - - ii (n, r?, r2, r) sin |^ — dr. 
Accordingly, the solution of the problem (10), (14) is 

p(ri,r2) = (ri,rg,r2,r^) sin 1^ ^ ^ / 2 Jo ('^i' ^i' ^2, r) cos |^ — dr, 

R (^1, T]*, r2, Tg) = Iq (^\/ (ri — r^){r2 — r2)^ is the modified Bessel function of the first kind of 
a zero order, having the following properties 

/„(o) = i, /;(o) = o. = ^ = 

Let us solve one practical problem to show the effectivity of using conservation laws instead 
of linearization by hodograph transformation. 

In work [14] the Cauchy problem for the cavity in an infinite medium loaded by a constant 
shear stress in addition to a uniform pressure (the so called Mikhlin problem) was solved under 
the condition Ji 7^ 0. Using conservation laws permits to forget about this condition. 

Let us consider an example, namely, let the contour be given by the 27r-periodic curve 

x = — rcott, y = —r, ^£[7 — vr, — 7), 

x = — rcott, y = r, tG(7, vr — 7), 



L 



vrt vrt 
^'cos— , y = rsm— , ^£[-7,7], 

27 27 

vrt TTt 
^^^^^^ \~^^ y = ?'sm— ^, t G (7 - 7r,7r - 7 , 

2(7r — 7) 2[TT — 7) 
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Figure 3. Characteristic field for tlie cavity. 



where a is the distance from (0,0) to the center of curve parts with radio r, 7 = arctan(r/a). 
This contour is loaded by the normal and tangent stresses 

0"n = -P, Tn = 0, 

where n is the normal to the contour. Let us put p = k = 1/2 and define the angle between the 
tangent to the contour and x-ax.es as N{t) = arctan |f, then the initial conditions take the form 

'iV(t) -7r/4 + ^/2, tG(0,7r), 

iV(t) -7r/4 + 37r/2, tG(7r,27r), 

-7r/4, t = 0, 

37r/4, t = vr. 

The solution of this problem is given in Fig. 3. Here a = 4, r = 3 and the first family of 
the characteristic curves is constructed, using conservation laws described above. Due to the 
symmetry of the contour, only the part with x > is shown. 

Let us note, that the stress field near the straight-line border of the cavity (regions PQM 
and P'Q'M') can not be constructed as a solution of a linearized system because Ji = in 
this domain. Moreover, the ordinary way to solve this problem is to consider the different 
parts separately. Firstly, one has to construct the field in the region PQM, assuming Ji = 0. 
Then one need to solve the Cauchy problem for the QQ' segment of the cavity. Having the 
solution of above problems, one should construct the field in the MQQ" region as a solution 
of the characteristic boundary problem, when the initial data are given along two characteristic 
curves MQ and QQ" . Using the conservation laws permits to avoid such a division of the 
boundary line and to solve the problem of Cauchy for the whole contour. 



4.2 Plane plasticity with Coulomb yield criterion 

The generalization of (21) is a system of the plane isotropic soil plasticity equations under 
Coulomb yield criterion, which in the case when the weight of the soil is neglected, has the 
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form [5, 26] 

— (1 + COS 2a cos 20) + — — cos 2a sin 20 
ox ay 

— — sin 20 — — cos 20 I , 
ox oy J 

— cos 2a sin 20 + ^r— (1 — cos 2a cos 20) 
ox oy 

— ^ cos 20 + sin 20 ) , 
ox oy J 

where 7r/2 — 2a is a constant angle of internal friction, a € (O, |), = + 7r/4 and k denotes 
the cohesion of the soil. If a = vr/4, then we have system (21). 

The Riemann invariants and the corresponding eigenvalues are the following ones (a 7^ 7r/4) 

a = u, 9 = V, Ai^2 = tan(0 it a), ri.2 = — - — ln(cj cot 2a + A;) it 0. 
If we take 

i?i = -n, R2 = -r2, ip{ri,r2) = y{Ri,R2), ?^2) = -R2), 

then equations for conservation laws (10) 

dip . f n-r2 \ dip dip fn-r2 \ dip 

coincide with linearized ones (8) 

dy (R1-R2 \ dx dy (R1-R2 \ dx 

tan a 7—— = 0, -— tan h a 77-— = 0. 



dRi \ 2 ) dRi ' dR2 \ 2 J dR 

This is a simple example, when one can need a change of independent variables ri, r2 to relate 
the conservation laws to the solution of a linearized system. 

Let us note, that introducing new functions $(ri,r2) and ^'(ri,r2) in such a way 

gcot 2a(ri+r2)/2 

ip = (-$ cos(0 + a) + ^' cos(0 - a)) , 

sin2a 

gCot2a{ri+r2)/2 

il) = (-$ sin(0 + a) + ^ sin(0 - a)) 

sm 2a 

one can reduce system (22) to the following one 

d^ ^ d^ ^ 

7^ + ^ . ^ =0, 7^ + ^ ■ ^ =0 
or2 2 sm 2a ori 2 sm 2a 

and finally to the telegraph equation of the form 

0. (23) 



dridr2 4sin^2a 
The corresponding initial conditions (11) look like 

= cos (!i^+„). 1^ 



(24) 



r2=r2 
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and for (14) we have 



^U=,o = -e-^(^?+-)sin(!i^ + a 



where 7 = ^ cot 2a. 

The solution of problem (23), (24) has a form 



dri 



0, 



r2=r§ 



(25) 



-^ir^^) cos 



+ a R{ri,r\,r2,r^) 



r2 



R{ri,ri,r2,T)e 



-7(r?+r) 



7 cos 



rV — r 



+ a 



sm 



rV — T 



+ a 



accordingly, the solution for (23), (25) is as follows 



sm 



^0 _ „0 
' 1 ' o 



+ a ) i?(ri,r5',r2,r^) 



ri — T 



^(ri,r°,r2,r)e-^(^?+^) 
where i? (n, r?, r2, r^) = Iq ( ~ ^i)(^2 - r^) 



7 sm ( -^—^ ^ ) ^ 2 



ri — T 



+ a 



4.3 Nonlinear hyperbolic heat equation 

As indicated in [8], the hyperbolic heat equation 
d fdU U\ d fxl dU' 



dt\dt ^ To) dX [mdX^^ °' 

where xo and tq are positive constants can be expressed in the form of the following quasilinear 
system 

du dv ^ dv Xo du ^ 
dx dy ' dx v? dy 

by introducing the potential function v{t,X) and setting 

U = Toe-^/^°u, y = X, x = e*/""". 
In this case 



ri2 = ue 



^"Z'^", Ai = -A2 



Xo 



n > 0, 



so ip = x(ri,r2), V = -yiri,r2)- 

Let us note, that in this case, equation (16) can be reduced to the following one 



d'^(p dip dip 



dxidx2 dxi dx2 



0, 



introducing new variables = e^^*. 

Explicitly, solution of problem (10), (11) looks like this 

2 3 1 a$ 

(p=—rlr^^, '4) = 2ri- h 

Xo ori dri 
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rir2 



4 h-o 

2 



/o(-ln2 3ln2- 



For problem (10), (14) function $ is the following one 



Xo 



1 1 i ^'i , I r2 
- m2 — ln2 — 

2 r? rO 



- In2 ^ ln2 — ] t idt 



$ = - 

(r0rir2) * 
4.4 Gas dynamics 

The one-dimensional isentropic flow of polytropic gas in Euler coordinates in case of plane 
symmetries, as is well known [21], is described by the following hyperbolic system 



ds , ^ .ds 

- + (a. + /3r)- = 0. 



dr , ^ ^dr 

- + (ar + /3.)-=0, 



where s = ri, r = r2 are the Riemann invariants, so 

Ai = an +/3r2, A2 = ar2 + /3ri, 

a = 1/2 + (7 — l)/4, /3 = 1/2 — (7 — l)/4, 7 = const 7^ ±1 is a parameter of polytrope. 
Explicitly, the solution of the problem (10), (11) has the form 



Xipi - X2P2 
Ai — A2 



Pi - P2 

Ai — A2 ' 



where 



P2{ri,r2) 



r2 - ri dp2 
K dri 



K 



7 + 1 
2(1-7)^ 



(rO - ri)^(r2 - r^yK+i) ^ l^i' ^2, ^2; 
K + 1 (t - r-?)^^ 



+ 



(r2-r?)(^+i) (t-n) 



F(r?,t;ri,r2)di, 



and F(r^, r^; n, r2) = 2-^1 (^-?^, + 1; 1; (^^_^o)(^o_^^) j is a hypergeometric function. 
Analogically, for problem (10), (14) function p2 undergoes a little modification 

, (M + arO)(rO -r0)2^+i n 



+ 



1 



(r2-r0)(^+i) 7.0 (t-n) 



(26) 



(r0-ri)^(r2-r0)(^+i) 

^ K F{rl t; n, r2) ( /3 + f^if^t + ar?) ) dt. 



In this case, equation (18) is satisfied, because h^^^ = /i^*^) = —a/3, so there is a relation be- 
tween the conservation laws and solutions of linearized system of the form ip = w{ri, r2)x{ri, r2), 
where function w is as follows 

w = {n- r2)(-+^)/(°-« = (n - r2)2/(^-i). 
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4.5 Loaded homogeneous semi-infinite elastic-plastic beam 

The process of propagation of plastic deformations in semi-infinite elastic-plastic beam, dyna- 
mically loaded on one end, in Lagrange coordinates is described by the following system [18] 

dv da dv 1 da 
^ dt dx^ dx pa?{a) dt^ 

where p = const is a density, v{t,x) is a velocity of medium particles, a tension a = a{e) 
is monotonically increasing convex function of the deformation e{x,t). Introducing function 
u{x,t), so that V = e = we come to nonlinear wave equation 

^-«^(-)^ = o> 

so a(e) is a speed of the longitudinal wave propagation in the beam. Note, that the above 
equation is of type (20). 

Eigenfunctions have the form Ai_2 = ^^(o") and Riemann invariants are the following ones 



ri 9 = v± 



" dai 
a{ai)' 



This is the case when Ai = — A2, then Lp = t{ri, r2), ip = —x{ri, r2). 
For example, if o = -y/a, then 



and system (10) can be reduced to Euler-Poisson-Darboux equation in the form 



dr\dr2 r\ — r2 \ dr\ dr2 

with the well known Riemann function. 

Moreover, this case is reduced to the system (26) with 

s = -ri, r = -r2, = -1/2, a = -1/4, /3 = 1/4. 

Solution of the problem (10), (11) has the form 

V' = 2 , V'= , Pi = /02 - 2(ri - r2) — , 

ri — r2 2 ori 

P2{n,r2) = ("-i^y F{rlrln,r2) - ^ ^ T ^^^F{rUri,r2)dt, 

\r\-r2J 2y^r^-r2JrO r^-t ' 

where F{r'(,r^;ri,r2) = 2-^1 (-1/2, 1/2; 1; (|^|I^O)(^^I^oj ) is a hypergeometric function. 
Function p2 for problem (10), (14) is the following one 



4 
3 



rr2 

/ Vn^F{r^^,t;n,r2)dt. 
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4.6 The Born— Infeld equation 



Let us consider one of the representations of Born-Infeld [4] nonlinear electrodynamics model 
in the form of linearly polarized plane wave equation 



1 + 



\ dx 



w 



. dw dw d"^ 



w 



dx dt dxdt 



1- 



y'dt 



w 



dx'^ 



0, 



which is a hyperbolic one for the solutions with l + ~ > ^- Introducing 



u 



dw 

dx ' 



dw 
~dt 



one can obtain the system in the form (4) 

du dv ^ dv f 2 — 1 du 2uv dv 
dt 1 + u'^ dx 1 + u'^ dx 



0. 



dt dx 

The Riemann invariants and eigenvalues are the following ones 

-uv + ^ ^ 



-uv 



ri 



and the corresponding system has the form (ri = r, r2 = 



Ai = r2, 



dr dr ^ 
dt dx 



ds f^-s _ Q 
dt dx 



-1, a = 0, /? = 1, i.e. the so-called Chaplygin 



The above system is a special case of (26) with 7 
gas. 

Note, that this case is an exceptional one, because Laplace invariants are zeros both for (15) 
and for (16), which take the form 



d'^x 
dridr2 



0, 



1 



d(p 



+ 



1 



dip 



d^v 

dridr2 ri — r2 dri ' ri — r2 dr2 



0. 



Solving (19), we obtained the relation of the components of conservation law if, tp with the 
solution of linearized system 



^{n,r2) 



x{ri,r2) 
ri - r2 



V'(ri,r2) 



2/(n,r-2) 



Let us note, that the above conservation laws are the particular case of ones calculated in [17]. 
Explicitly, the solution of problem (10), (11) has a form 



1 



r — s r — s 

and for the problem (10), (14) we have 
s , sr 



(27) 



(28) 



The final remark is that solution (12), (13), corresponding to (27) and (28) accordingly, is 
equivalent to the one provided in [2] under the condition for Jacobian Ji 7^ 0, and expressed 
in [17] in terms of Riemann invariants for the initial value problem: r\t=Q = r^{x), s\t=o = s^{x). 
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5 Conclusions 

Application of Riemann method to the system, describing conservation laws of quasilinear sys- 
tem, instead of linearized system, permits to solve the boundary value problems in the domain, 
where Jacobian of hodograph transformation is equal to zero. The proposed method permits 
the construction of characteristics, including both straight lines, corresponding to the simple 
wave solution (the constant solution too), and the case, when the initial curve is an envelope of 
a family of characteristic curves, which is useful in applications. 

In our paper is shown how the conservation laws (of the special form) of the considered 
system could be applied to construct the solution of a Cauchy problem, without separation of 
regions with a simple wave (or constant) solution, i.e. permit to find the point of characteristic 
without any restriction: is this point of straight line or not. As far we know it is a new point in 
the theory of considered form of quasilinear systems. 

The obtained theorem relates the Cauchy problem for quasilinear system with the solution 
of Cauchy problem for linear system of conservation laws. 

The relation between conservation laws of a quasilinear system and the corresponding linea- 
rized system, obtained by the hodograph transformation, gives nonsingular solutions. 

Recently, in [8] some non- homogeneous quasilinear systems, arising in various areas of physical 
interest, were reduced to the homogeneous ones by using the invariance to suitable Lie groups of 
point transformations. The well-known relation between the symmetries and conservation laws, 
we hope, opens the way to apply the proposed method for some non-homogeneous systems. 
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